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Abstract. The paper is devoted to homogenization of two-phase incompressible viscoelastic flows with 
disordered microstructure. We study two cases. In the first case, both phases are modeled as Kelvin- Voight 
viscoelastic materials. In the second case, one phase is a Kelvin- Voight material, and the other is a viscous 
Newtonian fluid. The microscale system contains the conservation of mass and balance of momentum 
equations. The inertial terms in the momentum equation incorporate the actual interface advected by the 
flow. In the constitutive equations, a frozen interface is employed. The interface geometry is arbitrary: we 
do not assume periodicity, statistical homogeneity or scale separation. The problem is homogenized using 
G-convergence and oscillating test functions. Since the microscale system is not parabolic, previously known 
constructions of the test functions do not work here. The test functions developed in the paper are non-local 
in time and satisfy divergence-free constraint exactly. The latter feature enables us to avoid working with 
pressure directly. We show that the effective medium is a single phase viscoelastic material that is not 
necessarily of Kelvin- Voight type. The effective constitutive equation contains a long memory viscoelastic 
term, as well as instantaneous elastic and viscous terms. 



1. Introduction 

Formulation of constitutive equations of multiphase materials under flow is a fundamental problem of 
continuum mechanics. Using mathematical homogenization theory [5], [15] . |25| to solve this problem is 
intuitively appealing, but not easy. The reasons for this can be summarized as follows: it is difficult to ho- 
mogenize evolution equations, non-linear equations, and equations involving general geometric distribution 
of the constituents. In particular, addressing the last difficulty is the necessary first step in developing a 
complete homogenization theory for moving interface problems. 

Consider a composite material with two constituents which we call phases. During flow, the interface 
between the phases is advected by the flow velocity. Therefore, the interface motion is coupled to the flow 
dynamics. A priori, one cannot expect that a geometry that is random homogeneous at the initial time 
remains random homogeneous at future times. Scale separation also cannot be expected to hold for all times 
in the interval of interest. Therefore, homogenization techniques that require a specific type of geometry, e.g. 
two-scale convergence ([ll[21]) and ergodic theorems, cannot be used in general. This leaves G-convergence 
[201 1^31 [571 [551 [311 [32] and, whenever variational formulation is available, F-convergence ([51[H]). The problem 
studied in this paper is not variational, so G-convergence is chosen as the main technical tool. 

G-convergence is a general notion of functional analysis and operator theory. In fact, G-convergence 
of an operator sequence can be identified with the Painleve-Kuratowski set convergence [THl [M] of the 
corresponding operator graphs F"^ (see e.g. the definition of G-convergence in [23]). Set convergence is 
sequentially compact provided the topology of F*^ has a countable base. Therefore, existence of an abstract 
G-limit operator A is easily obtained for A"^ that may be non-linear, non-local and multi-valued. Once 
existence of A is established, one needs to describe the structure of A, which is a problem of characterization. 
To solve this problem we use the method of oscillating test functions [51 [TUl HO] [53] . Our point of view is 
somewhat different from the standard one. To explain this further, consider a sequence of problems 

A'u' = f, 

where A'^ : X ^ X* are linear operators, X is a separable Banach space and X* is its topological dual. 
Suppose that A"^ G-converges to an operator A. This means that u'^ weakly in X and Au = f for each 
/ €E X*. Instead of choosing a sequence of oscillating test functions, we prescribe a sequence of corrector 



2000 Mathematics Subject Classification. 35B27, 35Q30, 35Q35, 74D05, 74E30, 74Q10, 76T99. 
^not necessarily random homogeneous 



1 



operators N"^ : X ^ X and then define the test functions by 

where w is an arbitrary fixed test function from a dense subset of X. 

It turns out that the characterization problem can be solved if N'^ satisfy 

(1.1) N^w w weakly in X, 

(1.2) (^^)"^ N^'w ^ A^w strongly in X*, 

possibly along a subsequence. Here, {A'^ j^ is the formal adjoint of A'^ . The operator aJ must be described 
as explicitly as possible. 
Indeed, if (fLTj) holds then 

(1.3) {A'u',N'w) = if^N'w) ^ {f,w) = {Au,w), 
and if (|1.2p is true, then 

(1.4) {A''u',N'w) = {u'',{A^)'^ N^w) (u,A^w) = (Au,w), 
Since w is any function in a dense subset of X, (|1.3p . p.4p imply Au = Au. 

In practice, finding a corrector operator is a difhcult task, and at present there are no general recipes 
for doing it. Known methods of G-limit characterization are specific to a certain class of problems: elliptic 
equations in divergence form ^20J , elliptic and parabolic problems with coercivity [3TJ [32] , linear elasticity 
[22], operator equations with monotone [7], and pseudo-monotone [53] operators. 

In most of these cases, the G-limit problem has the same general structure as the e-problems. For example, 
mixing two linear elastic materials produces an effective material that is also linear elastic. However, in 
mechanics of multi-phase flows, the phases may have different physical properties and the general nature of 
the effective equations is often unknown. Even in the relatively simple case of dense suspensions of rigid 
particles in a Stokes fluid, effective equations are subject to debate ([H[T31[1^). Therefore, the first goal of 
homogenization would be to determine the general structure of the effective equations: one-phase or two- 
phase, simple or not, type of memory dependence, presence of additional state variables etc. G-convergence 
is particularly well suited for answering questions of this kind. 

In this paper, we use G-convergence to study homogenization of two-phase incompressible viscoelastic 
flows. The phases are modeled as Kelvin- Voight materials: the elastic stresses satisfy Hook's law written 
in the spatial (Eulerian) formulation, and the viscous stresses obey Newton's law. Despite the fact that 
constitutive equations for the phases look similar, this model can describe the mixture of two materials, 
one of which is fluid- like (small elasticity), and the other is solid-like (small viscosity). We also study the 
fluid-structure interaction problem, where one phase is a Kelvin- Voight material, and the other is a viscous 
Newtonian fluid. At the initial moment of time the two phases are flnely mixed, and both phases occupy 
connected domains. The latter assumption excludes particle flows, since particles may collide, and that 
requires more complicated interface conditions. The density of each phase is constant, but the initial density 
of the mixture is highly oscillatory. This means that the sequence of initial densities does not converge 
strongly. No further assumptions are made about the initial interface geometry. In particular, we do not 
assume periodicity, statistical homogeneity or scale separation. 

The system of equations contains the mass balance and the momentum balance equations. The mass 
balance is needed because the initial density of the mixture is not constant, and the interface geometry 
changes in time. We also make the following choice regarding modeling of the interface. The mass balance 
equation and the convective terms in the momentum equation incorporate a moving interface advected by 
the flow velocity. The constitutive equations for the stress employ a frozen interface, that is the interface 
that existed at the initial time. This choice makes sense both physically and mathematically. 

From the point of view of physics, freezing the interface in constitutive equations is the only option 
compatible with the Hook's law. In linear elasticity, deformations are assumed to be small, and thus spatial 
and referential descriptions are identifled. In the referential description, the interface is always fixed. The 
reason for this is simple: a material particle that belongs initially to phase one remains in that phase at 
future times. Therefore, in the framework of linear elasticity, the interface must be fixed in both descriptions. 
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Otherwise elastic forces may lead to non-physical energy dissipation. More details illustrating this point are 
given in Appendix A. 

This may seem like a reason to use fixed interface in all terms and work with the formally linearized 
equations of motion. This is not satisfactory for viscoelastic composites, because convective effects may be 
significant even when the interface is nearly stationary. Indeed, when one phase is fluid-like, it may flow at 
large Reynolds number in the cavities of a matrix made from the second, solid-like phase. In that case the 
deformations inside the solid-like phase are small, as are the deformations of the interface. However, the 
deformations in the fluid-like phase may be large, and convective effects inside that phase cannot be neglected. 
Formal calculations, presented in Appendix [XI show that presence of convective terms is incompatible with 
the frozen interface. It is also impossible to meet the natural interface conditions unless convective terms 
are present in both phases, and not just in a fluid-like phase. 

To summarize, the microscale model in the paper partially accounts for the shape variation of the interface 
and satisfies all of the above natural requirements. In addition, these equations provide a good benchmark 
model for studying more complicated problems with moving interface such as flows of two immiscible fluids 
with surface tension, flows of nonlinear viscoelastic materials, fluid-particle flows etc. 

Mathematical study of the above microscale problem makes sense because little is known about G-limits 
of non-parabolic systems of fluid mechanics. Since the sequences of initial data do not converge strongly, the 
known compactness results for weak solutions such as Theorem 2. 4. in [18j and Theorem 5.1 in [19] cannot 
be expected to hold. In a sense, the problems of the type studied here are more difficult than existence 
questions. Because of the lack of compactness, the structure of G-limits for generic continuum mechanics 
systems may be different from the structure of e-problems, but this is not known. We show that in the present 
case, a new long memory term appears in the effective constitutive equations. Such results were previously 
obtained for linear equations with fixed interface. Formally this was done for locally periodic geometry in 
[25] and [6], and rigorously for periodic geometry in [12] using two-scale convergence. In [13], oscillating test 
functions were used for non-periodic scale separated geometries (some details of the proofs in that paper 
were omitted). In [12] and [13], the oscillating test functions did not satisfy divergence- free constraint. 
Using arbitrary test functions in incompressible problems makes it necessary to work with pressure directly. 
This is a serious obstacle, since good estimates on the pressure are not available even for density dependent 
Navier Stokes equations [18] . In this paper, we incorporate physically meaningful convective terms, allow for 
interfacial motions, and propose a construction of divergence-free oscillating test functions. The construction 
of the corrector operators involves certain auxiliary functions satisfying auxiliary problems. The right hand 
sides in auxiliary problems are chosen to satisfy condition (jl.ip on the corrector operators iV^ . The treatment 
here is inspired by the "condition N" in the papers [3I1I32] of Zhikov, Kozlov and Oleinik (see also the book 
[22j for an application to linear elasticity). It seems that many ideas in these works can be extended to 
non-parabolic evolution problems, including problems with moving interface. 

The main result of the paper are Theorem jSHj in Sect. |3l and Theorem lT.ll in Sect. |7| There we show that 
the effective system of equations describes a single-phase incompressible viscoelastic material. The effective 
system contains equations of mass and momentum balance. The effective constitutive equations contain 
a linearly elastic term, a linearly viscous term and a viscoelastic term that models a long term memory 
dependence of the effective material. 

The paper is organized as follows. Section 1 is introductory. Microscale problem formulation and prop- 
erties of finite energy weak solutions are given in Section 2. An outline of the existence proof for each fixed 
e > is presented in Appendix [Bj With some minor changes due to the presence of the elastic stress, we 
follow closely the existence proof for incompressible density-dependent Navier-Stokes equations in [18], Sect. 
2.3, 2.4. The global weak solutions of Leray type obtained in this way satisfy the energy inequality. Section 
3-6 contain a detailed study of the the case of two Kevin- Voight materials. In Section 3, we formulate 
the main theorem and provide an outline of the proof. Section 4 is devoted to constructing the corrector 
operators and oscillating test functions. Section 5 deals with passage to the limit in the inertial terms of the 
momentum balance equation. In Section 6 we obtain the effective constitutive equations and combine all 
the results to finish the proof of the main theorem. Finally, in Section 7 we indicate the changes necessary 
to treat the fluid-structure interaction case, and state the main theorem for this case. 
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2. Micro-scale problem. 



2.1. Equations of balance and constitutive equations. 

2.1.1. Choice of a model. We consider two-phase materials in which at least one of the phases resist shearing, 
and the material stress tensor can be written as a sum of an elastic (conservative) and dissipative stresses. 
To avoid mathematical difficulties in dealing with nonlinear elasticity we limit our investigation to flows 
for which the Hook's law of linear elasticity is an appropriate model of the elastic stress. We further 
suppose that deformation of the interface is small. In this case, the equations of motion are often formally 
linearized, assuming that the density is constant, and spatial and referential description are identified. As 
a consequence, the interface (always fixed in the referential description) is fixed in the spatial description. 
This is unsatisfactory when the physical properties of the phases are different., e.g. one phase is solid-like 
and the other one fluid-like. In this case, large deformations of the fluid-like phase may occur even when 
the deformations in the solid-like phase are small. Consequently, the contribution of the inertial terms to 
the overall momentum balance cannot be neglected. In addition, a correct model should respect physically 
realistic interface conditions (continuity of velocity and equality of tractions) as well as physically correct 
(formal) energy balance. These requirements make it necessary to consider inertia in both phases, with 
the corresponding terms that include a moving interface and densities governed by the mass conservation 
equations. By contrast, Hook's law for the elastic stress may lead to the non-physical dissipation of the 
elastic energy unless the interface in the constitutive equations is frozen. In Appendix A we present some 
(formal) calculations illustrating the above points. 

In this paper, we work with the micro-scale equations that employ a moving interface in the inertial terms 
and a frozen interface in the constitutive equations for the stress. This model can be viewed as a transition 
between a completely linear model with constant densities and static interface, and a fully nonlinear nonlinear 
that involves nonlinear constitutive equations for the elastic stress. 

2.1.2. Micro-scale equations. Let p'^,v' denote, respectively, the density and velocity of the composite. We 
also define 



(2.1) 



Jo 



Interface evolution equation. Evolution of V'^{t), W'^{t) can be described by the interface evolution 
equation. Let d'^{t,x) denote the characteristic function of V^{t). In the referential description, is fixed. 
Therefore, the material derivative of 0^ is zero. In the spatial description we have 

(2.2) dt9' -Ve' ^0 in U. 
Equation (j2.2p is supplemented with the initial condition 

(2.3) e^O,x) ^ e^g{x) inU. 

Mass conservation equation. Next, we state the mass conservation equation. The composite density 
satisfies 

(2.4) dtp' + divip^v') = 0, in U 
with the initial condition 

(2.5) p'{0, x) = p,0'o{x) + p2il - eiix)), xeU, 

where pi, p2 are the densities of the respective phases. These densities are assumed to be constant and 

bounded below by a positive constant. 

Incompressibility: 

(2.6) div v' = 0. 
Momentum balance for the composite: 

(2.7) dtv' + div{p'v' ® v') - div {T\ - Pf / + T| - P|/) 0. 
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Here — P^I is the stress tensor in the phase s, is the pressure, and / denotes the unit tensor. The 
initial conditions for (|2.7[) are 

(2.8) v^-{Q,x)^Vo{x), 

where v^) does not depend on e. In addition, (I2.1|l implies that u^{{),x) = 0. 
On the boundary dU, the condition 

(2.9) v^t,x)^0 
is imposed. 

Constitutive equations. As explained above, static interface seems to be a natural choice that is relatively 
easy to handle (unlike combining referential formulation for the elastic stress with the spatial formulation 
for the viscous stress), and compatible with the spatial form of the Hook's law for the elastic part of the 
stress. We therefore define 

(2.10) ^Tl+T^2- P"!^ where P' = Pf + P|, 

T\{t, x) = e'^{x) {A^e{u') + B^e{v')) , 

n{t, x) = {l- ei{x)) {A^e{u-) + B\{v^)) , 

and s = 1,2 are constant material tensors. In (|2.10p . e i(V + V^) is the symmetric part of the 

gradient. We assume that both phases are isotropic. In that case 

(2.11) Al^^i = fi^S.kSji, P,f,fe; = ,^'6,k5,i, 5 = 1,2, 

where fi" are the elastic moduli and i/" are the viscosities of the phases. All these constants are supposed to 
be positive. 

We also assume that the tensors A'^ , satisfy 

(2.12) ai^ • I < • ^ < -t Pi^-^<B'^-^< M ■ 
for each | G R^^^, with > 0, /3i > independent of e. 

Remark. Together, (|2.10p . (|2.4p . (|2.7p form a closed system, and thus the interface evolution equation (|2.2p 
can be dropped. This fact is important for compressible flows for which the mass conservation equation is 
stable with respect to weak convergence, while equation (|2.2p is not. For incompressible flows considered in 
this paper, the interface evolution equation has the same structure dt + div(i;'^.) as the mass conservation. 
Moreover, if the initial densities pi,P2 are constant, then the densities of the phases remain constant during 
the motion. In this case, the interface evolution equation and the mass conservation equation are essentially 
equivalent. 

Interface conditions. There are two interface conditions: the first is continuity of v across the inter- 
face (which is the actual moving interface governed by ()2.2p ). and the second is the equality of tractions 
(Tj — Pgl) Vs on the frozen interface. Here Ug denotes the exterior (to the phase s) unit normal to the 
frozen interface. 

2.1.3. Weak formulation of the micro-scale problem. In this section we provide the weak formulation of the 
problem to be homogenized. It consists of the mass conservation and momentum balance equations. 
Mass conservation. 

(2.13) / p^{Q,x)(t,{Q,x)dx- / / p'dtcfxixdt- / / p'v' ■ W(l>dxdt ^ 0, 

Ju J It Ju J It Ju 

where p^(0, a;) is given by (|2.5p . and 

/T = [o,r). 

Equation (|2.13p is supposed to hold for each smooth test function 0, equal to zero on dU and vanishing for 
t > T. 

Momentum balance. 

(2.14) -/ p''{Q,x)va-ij}dx- i / p^v^ ■ dt.ijidxdt - I i p^v" ® -Vijidxdt 

Ju J It Ju J It Ju 

+ I ( {T\+ Tl) ■ e{ifi)dxdt = 0. 

J It Ju 
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Equation (|2.14p holds for all smooth test-functions if), such that div -0 = 0, "0 equal to zero on dU, ip{t, a;) = 
when t >T. The dependence of on vf and 9q is given by (|2.10p . 

Remark. It is important to note that, because of the condition v'^ — on dU, the identity (|2.14p also holds 
for test functions ip with the condition ip — replaced by a less restrictive ip ■ v = on dU, v is the exterior 
unit normal to dU . This fact will be used in Section [3] to construct oscillating test functions. 

2.2. Finite energy weak solutions and bounds. We suppose that the initial conditions satisfy 

(2.15) < Ci < p'^(0,a;) < C2, 

(2.16) v^^HliU) 

with Ci, C2 independent of e. The system (|2.13l) . (I2.14[) closely resembles the system of density-dependent 
Navier-Stokes equations with density-dependent viscosity studied in [18 , ch. 2. The only difference is the 
presence of the strain-dependent terms of the type Ae{u'^) in the constitutive equations. When the viscosity 
does not vanish, as is the case here, the existence and overall properties of the weak solutions are determined 
by the viscosity, and not elasticity, of the medium. The proof of existence for each fixed £ > is outlined in 
Appendix iBl It yields, for each e G {£k}kLn existence of the finite energy weak solutions of (|2.13p . (|2.14p 
with the same properties as in [18j . theorem 2.1, namely 

(2.17) e L°^{It X [/) n C(/t, L^\U)), for aU 1 < p < 00, 

(2.18) €L\It,HI[U)), 

satisfying the energy inequalities 

(2.19) \ f p^v^'dx{t)+ f [A'd^,ix) + A'il-d^,ix))]\eiun\^dxit) 

^ JU JU 

+ [ [ [B^0'o{x) + B^{l~e'oix))]\e{v')\^dxdt< 

J It JU 

\ J^p^\v^\^dx{0) + J^[A%{x)+A'{l^eUx))] \eiu^)\^dxiO). 

Let us list the implications of the above estimates. First, renormalizing solutions of the mass conservation 
equations, as in [18j . sect 2.3, we obtain 

(2.20) \\P' \\l-(It>cU)<C2, 

Next we note that (|2.ip implies u'^{0,x) = 0, and by (|2.15|) . (|2.16p . the other initial conditions are bounded 
independent of e. This implies that the left hand side of (|2.19p is bounded independent of e, so that 

(2.21) \\T' \\L2iI^y,U)<C, 

with C independent of e. Combining (j2.2ip with the first Korn inequality for functions with zero trace on 
the boundary (see, e.g. [52], th. 2.1), and then with Poincare inequality, we deduce 

(2.22) II«MIl^(/.,h.5(c/))<C, 

(2.23) II \\L-ii^,Hiiu))< C. 
with C independent of e. Then it follows that 

(2.24) II p-\v-\^ \\L-(iTX-iU))< 
which, together with (j2.20p yields 

(2.25) \\p'v'\\LmT,L-{U))<C, 
with C independent of e. 



(2.26) 


P"(0,2 


(2.27) 




(2.28) 


— > 


(2.29) 




(2.30) 




(2.31) 


p'v' - 


(2.32) 





The above uniform bounds aUow us to extract a subsequence, still denoted by {sk]'k=i such that 

e) Pq, weak — * in L°°{U), 

p, weak — * in L°°{It x U), 

V, weak in L^{It,HI{U)), 

T, weak in L^{It, L^iU)), 

u, weak in L'^{It, H^iU)), 

^z, weakin {It, L^iU)), 

P, weak in {It , {U)) , 

The last comment concerns dtp^ and dt{p'^v'^). For technical reasons we need bounds on these sequences 
in the space L^{It ,W~'^^^{U)), where m > may be large. Such bounds can be also deduced from (|2.19p 
and (ITTa . (I2l4l) (see [19]). 

3. Main theorem and outline of the proof 
The main theorem is as follows. 
Theorem 3.1. The limits p,U, u satisfy 

(3.1) divtJ=0, 

and the integral identities 

(3.2) / po(/)(0, a;)da; — / / 'pdt4>dxdt — / / pv-Vi^dxdt — O, 

Ju J It Ju J It Ju 

(3.3) — / PqVq ■ xpdx — / / pv ■ dttpdxdt — / / 'pv^v-Vij^dxdt 

Ju J It Ju J It Ju 



I IT- e{ip)dxdt = 0. 
J It Ju 



I It Ju 

for all smooth test functions <f), tp, such that div ■0 = 0, and (f>,ip are equal to zero on dU and vanish for 
t > T. 

Moreover, there exist the effective tensors A e L'^{U),B e L'^{U) and C G L^{It x U) such that the 
effective deviatoric stress T satisfies 

(3.4) T = Ae(u) + Be(tT) + j C{t ^ T)e{v){T)dT. 

Jo 

Remark. Equations (|7.1iP " (|7.14p essentially mean that the effective equations for p, tJ, u, P are 

dtp + div(p v) = 0, 
div U = 0, 

(3.5) dt{pv) +div(p tJ(8) tJ) - div T + VP = 0, 
with initial and boundary conditions 

v{0, x) = vq, 'p{0, x) — 'Pq, tJ = on dU, 

and T given by ((711)) . 

The result follows from a number of propositions and theorems. Here, the outline of the proof is presented 
for the reader's convenience. 

1. Construct the corrector operators N'^w — . 

are defined by (I4.1|) - (l4.3p . The proposed construction is non-local in t and satisfies the divergence-free 
constraint. 

2. Auxiliary problems. Equation (|4.ip includes two types of auxiliary functions: n*"^'^ and mP'^''^ . An 
expression of the time derivative of w'^ additionally contains the final value mJ^'^ — m'P'^'^{T). These three 
families of functions satisfy auxiliary problems. The right hand sides in the auxiliary problems are chosen so 
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that the sequences of solutions converge to zero weakly in the appropriate Sobolev type spaces. The choice 
of the right hand sides involves abstract G-limit operators corresponding to each of the auxiliary problems. 
An analysis of auxiliary problems is presented in Sections 14. 11 14.21 Propositions 14. 51 14.81 are used repeatedly 
in the remainder of the proof. 

3. Convergence of w'^ . Using estimates from step 2, show that 

dtw' dtw weakly in L^{It, H^{U)). 

This is done in Proposition [5?2l This convergence corresponds to the condition (|l.ip on N^. 

4. Convergence of the inertial terms. Proposition 15.21 combined with the Lemma 15.31 (Lemma 5.1 from [18, ) 
implies that 

/ ^"(0, x)vq -w^+f f p^t!^ Wf+f I p'^v" ® ■ Vu;^ 
Ju Jo Ju Jo Ju 

converges to 

/ p{0,x)vq ■ w + / pvwt+ / pviSiv-Vw. 
Ju Jo Ju Jo Ju 

This is done in Proposition [HUl 

5. Compensated compactness of the stress. Convergence of the inertial terms implies 

r [ T^-e{w-)^ f I T-e{w). 
Jo Ju Jo Ju 

This is shown in Proposition l6.2l 

6. Effective stress. Characterization of T is obtained in Theorem 16. II This step corresponds to establishing 
condition (|1.2p on N'^ . The main idea is to write 

(3.6) / [ ■eiw') = - {u', div (A' -B'dt)w'). 

Jo Ju 

The expression for (u*^, div {A^ — B'^dt) w^) = (it^, div (A^ — B^dt) N^w) contains a number of terms de- 
pending on ■nP''-'^ ,m^'^ , mP'^ '^ and the corresponding pressures. In some of these terms we can pass to the 
limit using Lemma 5.1 from [TH], since it'^ has Sobolev regularity. In other terms this is not possible, but 
these terms vanish by design of the auxiliary problems. The effective tensors A, B and C are obtained as 
weak limits of the three fluxes that appear in the auxiliary problems for, respectively, •nP''''^, ml^'^ and m^'^''^. 

7. Mass conservation. Mass conservation equation is weakly stable. This is a well known fact (see |18|V 

4. Corrector operators and oscillating test functions 
We look for corrector operators of the form 

(4.1) N^w = w%t,x) = w{t,x)+nP'i^^{x)e{w)pq 

+ mP"-' {t-T + T, x)eiwt)pq{T, x)dT + V?i^ 

(Summation over p,q € {1,2,3} is assumed). Here, w G C^{It x J7), div u; = is an arbitrary test function, 
npq,e g Hq{U), mP'^'^ e L^{It, Hq{U)) are to be specified. So far we require that 

(4.2) div nP«'= = 0, div mPi'^ = 0. 

More conditions will be imposed below. 
The function (p" e L'^ (It , [U)) satisfies 

(4.3) Act)"" ^-nP'^'" ■\/e{w)pq- J mP"'' [t - t + T,x) ■ '^e{wt)pq{T,x)dT, 
V0= -1^ = on dU. 

The choice of the first three terms in (j4.ip is motivated by similar expressions used in periodic [12 and 
scale-separated |13| homogenization. The last term are need to enforce divergence-free constraint. This is 
necessary in order to avoid dealing with pressure in (|2.14p which is not L\^^(U) in general. 



Note also that Vcf)'^ ■ v is zero on dU for almost all t. This makes w'^ correctly defined test functions for 
(I2.14[) (see the Remark following that equation). Moreover, 

Lemma 4.1. The Junction defined by ^J^-^J^ satisfies 

div w'^ — 0. 

Proof. Taking divergence of (|4.ip and using (|4.2p we find 

div (^w{t, x) + nP'^'''eiw)pq + mP'^'%t -t + T, x)e{wt)pq{T, x)dT^ 

= nP'^^" ■ \/e{w)pq + mP^^'it -t + T,x)- \/e{wt)pq{T, x)dT, 
and the claim follows from the condition (|4.3p . 

■ 

4.1. Auxiliary problem for mP'^''^. In this subsection, p, q are fixed, so we drop them to simplify notations, 
and write instead of mP'^ '^ and so on. We look for that solve the auxiliary problem 

(4.4) - div {A^eirrf) - B^e{ml)) - VP| = /, div mf = 0, 
satisfying the condition 

(4.5) m^iT.x) ^m^'j.. 

The objective of this section is to show that the right hand side / can be chosen so that mf converges weakly 
to zero in an appropriate space. 

Let V(0 e C°°{It) satisfy ipiO) = and V^T) = 1. We use this function to reduce (gj]), (|43]) to a 
problem with a different right hand side and zero condition at t = T. Writing 

= rhf + ijjmfrp 

we deduce that rhf solves 

(4.6) - div (A^e(m'') - S^e(m^)) - VP^ = / + div [A'eiipmfp) - e{il)tmfp)) , div rh' = 0, 
with the inital condition 

(4.7) rh{T,x) = 0. 
Define the spaces 

(4.8) V = {ve L^{It, Hl{U)), div v = 0}, 

(4.9) W = {v eV -.vteV}, 

(4.10) Wt = {f e W : v{T) = 0}, 
The space V is equipped with a norm 

/ T \ 

(4.11) "^"^"(/ / e(^) • ■ 

This norm is induced by the norm 

(4.12) II V \\h},(u)= (^j^e{v)-e{v)dx^ ' . 

By Korn inequality ([12])j (|4.12p is a norm equivalent to the standard one. Also, Wt is dense in V. This 
can be proved in the same way as e.g. Thm. 2.1 in [17j. 

A weak solution rh,^ of (|4.13p , (|4.14p is an element of Wt satisfying 

(4.13) 



/ / (A^e(m^)-B^e(mD)•e(c/))da;d^ = (/-g^0> 
Jo JU 



for all </) G V. Here, 

= -div (A^e(V'm|,) - B^e(V'tm|.)) . 

Equation (|4.13p can be stated as 

(4.14) GW==/-s,^ 

with the operator : Wt — > V*. We consider it as an unbounded operator on V with the domain Wt- The 
corresponding bilinear form is defined as 

(4.15) {G^u,v) = [ I [A^ e{u) - e[ut)) ■ e{v)dxdt 

Jo Ju 

for each u G Wt, G V. Finally, we note that the adjoint operator G^'* with the domain 

Wo = {i; G W : v{{)) = 0} 

is defined by 



(4.16) 



(G^'*M,t;>= / I {A'e{u)+B'e{ut))-e{v)dxdt 
Jo Ju 



Proposition 4.2. (i) is strongly coercive: 

(4.17) {G'u,u)>ai\\u\\l 

where a\ is a constant from \2.1'2\) (and thus independent ofe); 
(a) G^ has a bounded inverse satisfying 

(4.18) ||(G-)-i/||^<^||/||^. 

ai 

for each / G V*. 

Proof. If u is sufficiently smooth then after integrating by parts in (|4.16p we would have 
(G^M, = y" ^'e{u) ■ e{u)dxdt B^e{u) ■ e{u)dx^ (0), 

where we took into account that u{T) = 0. However, for an arbitrary u G Wt the second term in the right 
hand side may not be well defined. To bypass this difficulty, observe that for almost all t G It, 



\(^j^B'e{u)-e{u)dx^ {t) 



is finite. For such t we have 



J (A^eiu) ■ e{u) - B^e(Mt) • e{u) dxdt 

= j A^e{u) ■ e{u)dxdt + \ (^J B^e{u) ■ e{u)dx^ (t) 

> A^e(tt) • e{u)dxdt > ai / e{u) ■ e{u)dxdt 

Jt Ju Jt Ju 

The last inequality follows from (|2.12p . Using absolute continuity in t of the first and last terms in the above 

inequality, we can pass to the limit t ^ 0^ and obtain 

(G'^tt, u) = lim I I {A^e{u) ■ e{u) - B^e{ut) ■ e{u) dxdt 
^^'^+Jt Ju 

> Hm ai / e{u) ■ e{u)dxdt = ai || m ||y, 
*^"+ Jt Ju 

which proves (i). 

(ii) follows from (i). This is known (see, e.g. [32], Lemma 1). We only sketch the proof for completeness. 
Since G^ is closed, passing to the limit in (|4.17p we obtain that the image of G"^ is closed in V*. If this image 
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does not contain all V*, then, because of density of Wt in V, there is g E V such that {C^u^g) = for all 
u e Wt- This yields g & Wo (domain of G"'*) and G^ *g = 0. Next we observe that G"'* satisfies (liTZ)) 
which yields g = and gives a contradiction. Thus G^ is onto. The estimate (|4.18p follows from (|4.17p . 

■ 

Remark, (ii) implies existence of the pressure e L^{It, LP'{U)). This follows using standard arguments 
from ,2H combined with the inclusion A^eljaP'') + B^e(mj'') e L"^ (It , (U)) . Moreover, P§ is bounded 
in (It , L"^ (U)) independent of e. Therefore, extracting a subsequence if necessary, we can assume that 
P| P3 weakly in L^{It, L'^{U)). 

Definition 4.3. We say that the sequence G-converges to an operator G : D{G) C Wt V*, if for 
each / e V* the sequence it"^ = (G')~ / converges to some u e D{G) weakly in Wt- In this case we define 
Gu = f. 

Proposition 4.4. The sequence G^ contains a G-convergent subsequence. The limiting operator G has the 

following properties: 

(^) 

(4.19) {Gu,u)>ai\\u\\l 

for each u G D{G); 
(ii) D{G)^Wt- 



Proof. Let us write G'^ ~ A' — B^dt where A'^ : V ^ V* are operators induced by the bilinear forms 

a^{u,v)= i i e{u) - e{v)dxdt, 
Jo Ju 

b%u,v)= I I B^eiu) - e{v)dxdt, 
Jo Ju 

respectively. EUipticity assumptions (|2.12[) imply that are coercive and bounded with coercivity 

constants and bounds independent of e. In particular coercivity of implies that there exists a bounded 
inverse (-B^) ^ defined on V*, satisfying || {B'^) ^ ||< where /3i is the lower bound from ()2.12|) . Therefore, 
if {A^ - B''dt)u'' = / then dtu^ = {B'T^ {A'u - f). This implies 

"2 II J, _c II ^ / "2 , II J. II 



dtu' llv<^ll (G^)"'/||v + ||/||v*<( 



O IIV" / J \ \ 1^ ' \\ J 111^ — \ o ^/IIJ III' ' 

Pi aiPi 
Thus, if m'^ ^ M weakly in Wt then 

(4.20) II w||w<G(ai,a2,/3i) || / ||v- 

Since V* is separable, we can use diagonal procedure to find a subsequence, non relabeled, such that = 
(G"^) ^ / converges weakly in Wt to u = G^^ f for all / in a dense subset of V*. Inequality (|4.20p implies 
that convergence also holds for all / £ V*, and the operator G~^ is bounded. 

Next, consider a sequence it*^ such that G^u'^ = f. Then, by the preceding, converges weakly to u, 
and Gu = / by definition of the G-limit. Since {G^u^,u^) — {f,u^) ~ {Gu,u^), we can pass to the limit 
and obtain 

lim(G''M^M'^) = (Gu,u). 

This, together with lower semicontinuity of the norm with respect to weak convergence, allows passage to 
the hmit in (|¥T71) which yields ^J^. 

To prove (ii), observe first that by (|4.20p . D{G) C Wt- To prove equality, we first show that G~^ is 
injective: G~^ f = implies / = 0. Arguing by contradictions, suppose that there is g € V*,g ^ such that 
G~^g — 0. Consider the sequence tt^ — (G^)~^g. By definition of G, it^ converges to zero weakly in Wt- 
Then by (|4l7l) . 

(g,«P-(G-tt^,ttP>ai im II?,. 



passing to the hmit e — we obtain that || \\v^ and thus Ug converges to zero strongly in V. Next, 
we use Ug ^ as the test function in the weak formulation of G^Ug = g. Integrating by parts (which can be 
justified as in the proof of Proposition 14. 2p and using coercivity of we obtain 

and thus Ug ^ converges to zero strongly in V. Next, using uniform boundedness of , B^ , we write 

II 9 \\v* = \\ A'ul - B'ul t ||v*< 02 II ul llv +/32 || ||v 

where a2, P2 are constants from (|2.12p . Passing to the limit in the above we deduce g — 0, which contradicts 
assumption g 7^ 0. Thus is injective. 

Next we show that D{G) (equivalently, the range of G~^) is dense in V. If this were false, there would be 
h ^ 0, h E V* such that {h, G~^f) ~ for all / £ V*. Let Uh = G~^h. Choosing f — h, we obtain using 
(I4l9ll : 

= {h,G^^h) = {Guh,Uh) > ai \\ Uh ||v • 

Therefore, Uh = 0. Then h = by injectivity of G^^. This contradicts the assumption h ^ Q. 
Finally, observe that the norm of Wt induces a scalar product 

{u,v)wt ^ / / {e{ut) ■ e{vt) + e{u) ■ e{v))dxdt. 
Jo Ju 

In search of a contradiction, suppose that D{G) is a proper subset of Wt- Then there is U ^ 0,11 € Wt 
such that (U, G^^/)wt = for all / G V*. Expression {u,v)wt defines a bounded linear functional lu{v) 
on D{G) which by Hahn-Banach theorem can be extended to a bounded linear functional on V and this 
extension has same norm as lu{v). Therefore, (U, G~^/)wt — implies (L^, G"""^/) = for all / e V*. 
Density of the range of G~^ implies = 0. But then (it, m)wt = which contradicts the assumption 
u ^ 0. Thus (ii) is proved. 

■ 

Proposition 4.5. There exists f G V* such that the sequence of solutions of j^.^p , with this choice of 
the right hand side, contains a subsequence (not relabeled) satisfying 

(i) weakly in Wt, 

(ii) strongly in L"^ {It , L"^ {U)) . 

Proof. By Proposition 14.81 proved below in Sect. 14.21 we can assume that nij, converges to zero weakly in 
Hq{U). Then ipnij. — *■ and dtipnij. weakly in V. 

Since m'^ = mf + i/imy, (i) will be proved if we show that there is choice of / such that — ^ weakly 
in Wt- 

Consider (|4.14p . In view of (|2.12p and uniform bounds on m|,, the sequence is bounded in V*. 
Therefore, the sequence (G'^)~^g^ is bounded in Wt, and we can extract a subsequence that converges 
weakly to some q G Wt- By Proposition 14.41 (ii), q E D{G). Therefore, we can choose 

(4.21) f = Gq. 

Then 

= (G^r' if - gn = icn^' Gq - {g^'' g^- 

By Proposition 14.41 (G"^) ^ Gq q weakly in Wt, up to extraction of a subsequence. Hence, rhf 0. 
Thus (i) is proved. 

To prove (ii), observe first that m|, converges to zero strongly in L'^{U) and thus -tj^mfp — > and dtipm^ 
strongly in L'^ (It , L'^ (U)) . Therefore, to prove strong convergence of m'^ it is enough to prove strong 
convergence of rh^ . 

Next, note that rhj is bounded in V independent of e. Now strong convergence of rh^ is deduced from (i) 
and J. L. Lions' compactness theorem (see e.g. [29], thm. 2.1, ch.III). 

■ 
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4.2. Auxiliary problems for nP'' '^ ,m!^''^ . In this section, p,q are once again fixed, so we drop them to 
simplify notations, and write n'^ instead of n^"^'^ and so on. 

We seek , mfp satisfying, respectively, the auxiliary problems 

(4.22) - div {A' {P'^ + e(n^))) - VPf = /i, 

(4.23) - div {B"" (F' + e (n^)) - e (m|,)) - VP| = /a, 

with suitably chosen /i, /2- First, we find nf from (|4.22p . Then this nf should be plugged into (|4.23[) . and 
then m?p can be found. The goal, as before, to choose /i, /2 so that ,m?p would have subsequences that 
converge weakly to zero. Let 

V = {v eH^ : div V = 0}, 
equipped with the norm (|4.12p . Given f i E V* , £ V is a. weak solution of (|4.22p provided 

(4.24) / A'{IP^ + ein'))-e{w)dx = {f„w)v,v* 

Ju 

for all w E V. Weak solutions of (|4.23p are defined similarly. This identity can be written as an operator 
equation 

(4.25) A'n' = f,-gl, 
where 

= -div(A^F')ey^ 
and A'^ : V ^ V* is the operator induced by the bilinear form 

a^{u,v)— / A'^e{u) ■ e{v)dx. 
Ju 

Similarly, (I4.23P can be written as 

(4.26) B'm'j. = /2 - 
where the operator is induced by the form 



b'^{u,v)= / B^e{u) ■ e{v)dx, 
Ju 

g| = -div(B^(JP^ + e(n^))). 



and 

By (pi^ . operators satisfy 

(4.27) (A'u,u)v,v* >ai\\u\\l. 

(4.28) \\ A'u\\v*<a2\\u\\v 

(4.29) {B'u,u)v^v*>f^i\\u\\l 

(4.30) \\B'u\\v*<f32\\u\\v 

Lax-Milgram lemma implies existence of unique solutions of (|4.25p . (I4.26p . These solutions satisfy 

= {A^)-' if, - g\), II ||v< - || /i - gl \\v*, 

ai 

m|, = {B^y^ if 2 - gl), II m^T \\v< II /2 - 9l Ik* ■ 

Pi 

Remark. Existence of the pressures Pf E L^{U), P| E L'^{U) follows using standard arguments from ^S]. 
Moreover, these pressures are bounded in L'^(U) independent of e. Therefore, extracting a subsequence if 
necessary, we can assume that P? Pj,j ~ 1,2 weakly in L'^{U). 

Definition 4.6. The sequence of operators A^ : V ^ V* G-converges to an operator A if {A^) ^ j converges 
to some u E V weakly in V, for each f E V* . We also define / = Au. 

Proposition 4.7. The sequences A'^^B^ contain G-convergent subsequences. 

Proof. This is known |20j . thm. 2. 
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Proposition 4.8. There exists fi G V* (respectively f^) such that, up to extraction of a subsequence, 
(respectively rrij,) converge to zero weakly in V. 

Proof. Consider (|4.25p . Since gl bounded in V*, and is bounded independent of e, the sequence 

(A^) g\ is bounded in V . Thus we can extract a subsequence that converges weakly in V to some Ui G V . 
Choose 

(4.31) /i=A«i. 
Then 

By definition of A, the first term in the right hand side converges to Ui weakly in V , and so does the second. 
Hence ^ weakly in V . For (|4.26[) the procedure is the same. Up to extraction of a subsequence, 
[B^ ) ^ g2 weakly in V, and we choose 

(4.32) /2 = Bu2. 



5. Inertial terms in the momentum balance equation 

In the remainder of the paper, we assume that the oscillating test functions are defined as follows. 

Definition 5.1. Let w E C°°{It x U),c\\vw = be arbitrary, and define by (|4.ip . In (|4.ip . choose 
n'P'i-'^ ,171^'^ that solve, respectively (|4.22p . (|4.23p with the right hand sides chosen according to (|4.3ip . (|4.32p . 
Also, let m^''^ satisfy (|4.4p with the right hand side chosen according to (|4.2ip . 

Proposition 5.2. The sequence defined as above satisfies 

(5.1) in L2(/t, L'^iU)), -> Wt in L'^{It, L^{U)). 
Also, e L°°{It,H^{U)), and 

(5.2) II \\L^iI^.Hiiu))< II < \\L^(I^,HliU))< CT^'^ 
with C independent of e. 

Proof. First we need a formula for the time derivative of . After taking time derivative of (j4.ip . integrating 
by parts in the time convolution (which involves putting time differentiation on Wt instead of mJ"^'^) and 
using w{T) = we obtain 

(5.3) wl = Wt+'nP''-'^{x)e{wt)-pq + 




mPi^^it -T + T)e{wu)pq{T)dT + V0f . 



To prove strong convergence of tuf we need to prove that all terms in the right hand side of (j5.3p converge 
to zero strongly in L"^ {It , L"^ {U)) . 
Step 1. Show that 

nP''-'e{w)pq + mPi-'it ~t + T)e{wt)pq{T)dT ^ 0, 

and 

nP'''%x)e{wt)pg + ^ mP^'^t -r + T)e{wu)pq{T)dT ^ 

strongly in L^{It x U). 

By Propositions 14.81 nP'^''^ weakly in Hq{U)) and thus strongly in L?{U). By Proposition 14.51 
mP'^'^ — > strongly in L^{It, L'^{U)), and we conclude. 
Step 2. Show that -> and dtV(j)'' strongly in L^{It x U). 

First we estimate e(V0'^). Note that e{V4)^)ij = didjcj)'^ and write didj<j)'^ — didjE-k{A(l)'^) + didjK-k{A4>^), 
where E = c-^ is a fundamental solution of the Laplacian, and if is a harmonic function, which depends 
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only on U. To estimate didjE * (A(/)'^) we use Calderon-Zygmund inequality (see, e.g. [TP, Theorem 9.9) 
with p — 2 and obtain 

II d^djE^{A4)^) ||l2(/^xc/)<|| WmiT-^u) 
Since didjK is a smooth function, there exists a constant C{U) depending only on U such that 

II d,d,K^{A<P') \\mi^^u)< C{U) II A<^^ \\L-(ir>^u) 

Thus 

(5.4) II d,d,(p^ ||l^(/,xc/)< (1 + C^((7)) II A<^^ lU 
Combining (|5.4p with (|4.3p we find 

3 

(5.5) / / |e(Vr)pda;di<C(«,,C/) ^ (||nP^-||i.(^) + ||mf^'^||i.(,^,^) 

•' It JU p,q=l 

so e(V(/)^) converges to zero strongly in L^{It x U). This also implies that the components of the Hessian 
of (jf converge to zero strongly in LP' {It x U)- 

Finally, the standard a priori estimate for the Neumann problem A(/)^ — f^, f ^ satisfying Wcf)'^ -v = 

on the boundary, yields 

(5.6) J^^cl^' -Vcl^'dx <\\ r \\l2^u)W \\mu)< C \\ r Wmu) (^j^Vcj)' ■V4>'dx^ ' . 

Here, C is the constant in Poincare inequality. Poincare inequality applies after we impose the condition 
^Qjj 4>^dS = 0, standard for Neumann problem. Since given by the right hand side of (|4.3p . converges to 
zero strongly in L^(/t, i^(J7)), (|5.6p implies that Wcj)'^ converges to zero strongly in L'^ {It , L'^ (U)) . 
Differentiating (|4.3p in t and integrating by parts as in (|5.3p we find 

(5.7) A0? = • We{wt)pg - £ mP''^'{t -t + T)- We{wtt)pg{T)dT. 
Therefore, arguing as above we have 

(5.8) 



/ / \e{V(j}l)\^dxdt <C{w,U) V] fll" 



Pij^e ||2 I II rnPq,^ ||2 



which yields e(V(/)f) 0, and then V(/)f strongly in L^{It x C^)- 
Step 4- Prove 



Since w'^{t) = — /^^lof (r)dr, we obtain for almost all t ^ It 



II IIh^Hc/) W < [ II '^t llff^Hf) i^)dr < T'/^ II wl ||^,(^) (r)drj 



1/2 



and (lOl) follows. 



Next we will need the following lemma ([H], lemma 5.1). 

Lemma 5.3. Let g",h^^ converge weakly to g,h, respectively in L^^ (0, T; L^'^(rj)), L'^ (0, T; L'^^(rj)), where 
1 < Pl,P2 < oo, 

1 11 1 _^ 
Pi qi P2 92 

We assume in addition that dtg"" is bounded in L^{0,T;W^"^'^{n)) for some to > independent of n and 

II - /i"(- +Ct) \\Lii{o,T:Li2(n))'^ 

as 1^1 — > 0, uniformly in n. 

Then g^^h" converges to gh in the sense of distributions on fl x (0,T). 

This lemma can be used to obtain the effective mass conservation equation and to pass to the limit in the 
inertial terms in the momentum equation (|2.14p . 
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Proposition 5.4. Let w'^ be functions from Definition \5.1\ Then 

(5.9) lim / {Q,x)vq ■ dx — I p{0,x) vq ■ w dx, 

(5.10) lim / p'^v^ ■ dtw^ dxdt = / pv-dtwdxdt 

J It X U J It X U 



an 



d 



(5.11) lim / p^ ® ■ V dxdt — / p v ® v -Vw dxdt. 

^^'^JiTy.u Jitxu 

Proof. By Lemma [5.31 p^v^ ^ p tT in the sense of distributions, and thus also weakly in L^{It x U). By 
(|5.ip . w^,dtw'^ converge to respectively w,dtw strongly in L'^{It x U). This permits passage to the limit 
in the products and yields (|5.10p . 
Since ~ w e C{It, L'^{U)), 

\\w^ -w\\l2[u){^) < -Wt\\L2(u) {r)dT 

Jo 

<T'^'( r \\wl-Wt\\mu){T)dT] . 





Noting that lof — *■ Wt strongly in L^{It x U), we obtain that (tt)'^ — ''^)(0) — ^ strongly in L?'{U). Since 
p^(0, x) converges weakly-* in L°°{U) to p(0, x), strong convergence of w^{Q, x) permits passage to the limit 
in the product ^"^(0, x)vq{x) ■ 10^(0, x) and yields (|5.9[) . 

Next, fix i £ {1, 2, 3}, pick a function 77 S C§°{It x J7), and insert the test function {w'^ — w)rj into the 
weak formulation of the mass balance equation. This yields 



(5.12) Jjjp^{0,x){w^ -w)r]{0,x)dx - / p' dt {{w' ~ w)ri) dxdt 

J It Ju 

- / / p^v^ ■ V {{w" - w)ri) dxdt = 0. 

J It -ju 

Strong convergence of (lo'^ — w){Q) to zero implies 

(5.13) lim / p''{0,x){w'' -w)ri{0,x)dx = 0. 

^^°Ju 



Next, note that 



(5.14) lim / / p'dt ((to" - w)ri) dxdt = 

2f 



because {wf — Wt) strongly in L {It x U), and p"^ is bounded in L°°{It x U) independent of e. Now 
from ([5T2I) . (jSTTS]) and ([STM]) we deduce 

(5.15) lim / / p''?;^ • V ((w^ - w)77)da;di = 0. 

J/t Ju 

Since p'^v'^ is bounded in L^{It x [/), and w'^ w strongly in L^{It x U), (|5.15p implies 

(5.16) lim / / -qp^v^ -Viw" - w) dxdt = 0. 

■^-♦o J It Ju 

Since 77 € C^{It x J7) is an arbitray test function, p'^i)^ • V (w^ — ii;) — > in I?'(/t x U). 

Next we claim that p'^v'^ ■ V (lo^ — w) is bounded in L-^{It, L^^^{U)) independent of e. This follows from 
Sobolev imbedding for v'^ and Holder inequality. Application of Holder inequality yields 

\p^vid,w^^\' dx^ (t) <\\ p- iii^(^) (t) (^J^ \vir^ ' it) (^j^ \d^wr'''Y {t). 
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Here s.q > 1 and - + = 1. Hence, 

71 — q q' ' 

2 



Therefore 

(5-17) II p^vldkw] ||L2(j^^i,(t;))<|| ||l-(/txc/)|| 9fcw| |lioo(i.,')|| v| \\l^It,l^-'{u)) 

We need to choose s, q so that (i) the right hand side of (|5.17p is finite; and (ii) G L^{It, (U)), where 
■i + ^ = 1. By Sobolev imbedding, s' < 6, and therefore 

(5.18) s>^. 

5 

By (Eai, II dkw'^ Ill5l(i.,') is finite if 

(5.19) 1 < sq' <2^=>1-- <s<2--, 

q q 

[sq' < 2 are aUowed because U is bounded). Also, by Sobolev imbedding 

1 6 

(5.20) 1 < sq < 6 - < s < -. 

q q 

The solution set of inequalities (|5.18p ~ (|5.20p is a non-empty, convex quadrilateral in the \/q — s-plane. For 
example, we can choose s = f and any q satisfying ^ < ^ < |. If, for example, 1/q — |, then sq — 3, 
sq' = I and (|5.17p becomes 

II p^V^-Ww^ ||^,(,^_^|(^^^< C II p^ \\L^iI^.U)\\ ILoo(,,,^f llL=(/.,L3(r/)) . 

Since the right hand side of (|5.17p is bounded independent of e, the claim is proved. 

Together with ([Sle]) . this yields p^v" ■ V(u;^ - w) ^ weakly in L'^{It, L'^^^iU)). Therefore, the 
weak limit of p'^v'^ ■ Vw'^ is the same as the weak limit of p'^v'^ ■ Vw in L^[It,L^^^{U)). By Lemma [5.31 
p'^v^ ■ Vw fm ■ Ww. Thus 

(5.21) p^v" - Vw" -^-pv - Vw 

weakly in L'^{It, L^^^{U)). The bound on dt[p'^v'^ ■ Vw"^] in a negative Sobolev space follow from the 
corresponding bounds on p^v'^ and the fact that wf is bounded in L^{It, Hq{U)). 
Now application of Lemma [5.31 with — p'^v^ ■ Ww^, /i^ — v'^ yields 

(5.22) lim / p'^v^ ® ■ Vw^rj dxdt = pv(^v- Vwr/ dxdt. 

J It X U J It x U 

For each rj S C^{It x U). From Sobolev imbedding and bounds on Ww^, in the same way as (I5.17P was 
analyzed, choosing 1/q — 2/5, s — 6/5, we obtain 

(5.23) II p'v'(Sv' ■ Vw' ll^2(^^.i§ ([;))< C II p' \\l'^{ItxU)\\ Vu;^ \\L'^iiT,LHu))\\ v'(g,v' \\miT,L^u)) 

with C independent of e. Note that the right hand side is bounded independent of e. Therefore, p^v"^ 

■ Vw^ converges lo'p v ®v ■ Ww weakly in L'^ {It, L's {[/)). This implies convergence of integrals of 
p'^v^ (g) ■ Ww^ over subsets of It x U, and in particular (|5.1ip . 



17 



6. Effective deviatoric stress. Proof of the main theorem 

Theorem 6.1. There exist a subsequence, not relabeled, and effective material tensors A G L'^{U),B £ 
L'^{U) and C G £^(/t x U) such that for each w € C^{It x U) with div w = Q, w{T, x) = 0, 

/ T ■ e{w)dxdt = lim / T" ■ e{w'')dxdt 

= I I ('Aeiu) +^e{v) + [ C{t - t, ■)e{v){T, ■)) ■ e{w)dxdT 
J It Ju \ Jo J 

as e along this subsequence. 

Proof. The theorem foUows from several propositions. First, we prove that convergence of inertial terms 
impUes compensated compactness of stress. 

Proposition 6.2. Let be test functions from Definition \5.1\ Then 

Van i / • e{w^)dxdt = T ■ e{w)dxdt. 

J It Ju J It Ju 

Proof of the proposition. First, use ii; as a test function in (|2.14p and pass to the limit e — > 0. Repeated 
application of Lemma 15.31 in the inertial terms yields 



(6.1) 



/ p(0, x) Vq ■ w dx — p V ■ dtw dxdt 

Ju J It X U 

/ 'pv i^v ■ Vw dxdt + / T ■ e{w)dxdt = 

J It X U J It X U 



IItxU JItxU 

Then insert lo^ into (|2.14p and pass to the limit e ^ 0. By Proposition 15.41 the integrals corresponding to 
the inertial terms will converge to the corresponding integrals of the limiting functions p,v. This yields 



(6.2) 



- / p(0, x) Vq ■ w dx — p V ■ dtW dxdt 

Ju Jit/U 

I 'p V ®v -Vw dxdt + Van I ■ e{w^)dxdt = 

JitxU '^^^JitxU 



Comparison of (|6.ip and (|6.2p finishes the proof. 
Next, using symmetry of A'^ , we have 



[ ■ e{w^)dxdt ^ [ e(M'^) • [A^e(u;^) - B^'eiwl)] dxdt. 

J It X U J It x U 

Differentiation in (|4.ip yields 

(6.3) A'e{w')-B'e{wl) = .Fr^e(u;)p, + .Ff '=6(1^*)^? + -^f'^'l^ - ^ + r)e(«;t)p,(r)dT 

+A'{gl+gl)-B'{dtg\+dtgl) 
+A'e{V(j)') - B'e{V<Pl), 

Tl'^^'ix) = A' (F« + e (nP«-^)) , F« = ^ (e^ ® e, + e, ep) , 



where 



.Ff '"(a;) = B^ (JP* + e («''«'") - e (mP?'")) , 
Tl'^'^t, x) ^ A^e (mP«'") - B'e (mf ^") , 

(6.4) g\ = ]^ (nP«'-^ ® Ve(io)p, + (nf«'^ ® Ve(i«)p,)^) , 

(6.5) »2 = ^ ^ (mf9'^(t -T + T)® We{w)pq{T) + (mP«'^(t - r + T) ® Ve(t(;)p,(T))^) 



dr. 
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Next we show that the only terms in (|6.3p that contribute to the effective stress are the terms containing 
Proposition 6.3. Let he as is Definition ] 5. 1\ Then 

(6.6) lim / / e{u') ■ [A'{g\ + gl) - B'{dtg\ + dtg'^)] dxdt = 0, 

(6.7) lim / / eiu^) ■ [A^e(V^") - B"e(V(/)f )] dxdt = 0. 

J It Ju 

Proof of the proposition. Since e(M^) converges to e{u) weakly in L^{It x U), it is enough to show that all 
terms in brackets in (16. 6|) . ()6.7p converge to zero strongly in L?{It x U). 
Step 1. Prove fO)) . 

By Proposition 14.81 n^'^^^ converges to zero strongly in L?{It x U) Therefore, gi,dtgi in (|7.3p converge 
to zero strongly in L^{It x U). By Proposition 14.51 mP'''^ converge to zero strongly in L^{It x U). Hence, 
g| in (j7.4p converges to zero strongly in L^{It x C/). Next, differentiate g| in i and integrate by parts in the 
time convolution exactly as in the proof of Proposition [5T2l Then 

dtgl (mf«'=-(t -T + T)® We{wu)pq{T) + (mP«'^(i - r + T) ® Ve(-«;„)p« (t))^) dr. 

which converges to zero strongly in L^{It x U). Next, since A'^,B'^ are bounded pointwise independent of 
£, we deduce that 

A'{gl+gl)-B%dtgl+dtgl) 

converges to zero strongly in L^{It x U). 
Step 2. Prove [6l^ . 

From we have e(V(/)^) 0, and by jEH]), e(V(/)f) ~* strongly ^^(/t x U). Hence, A^e(V(?!)^), 

-B^e(V0f ) also converge to zero strongly in L^{It x U). 



By Proposition!; 

/ T ■ e(w)dxdt = lim / / • e{w'')dxdt = lim /(«", u;^). 



where 



(6.8) I{u',w') =11 e{u') ■ ( J^f ■^e(-u;)p, + J^'^^' e{wt)pq + f J^'^^'it - r + T)e{wt)pq{T)dT\ dxdt 
J It Ju \ Jt J 

Proposition 6.4. There exist the effective tensors A E L'^{U),B E L'^iU) and C E L^{It x U) such that, 
up to extraction of a subsequence, 



(6.9) lim£^o^(M^,t"^) = -{u,dW (^Tl'^e{w)pq^ +diy (^Tl'^eiwt)^ 

-(M%div Tl'^{t-T + T)e{wt)pg{T)dTdxdtj) 

^ f f (^e{u) +Be(ut) + [ C{t - t, ■)e{ut)iT, ■)dT\ ■ e{w)dxdt, 

J It Ju \ Jq ) 

Proof of the proposition. First, we note that J^j"''^, Pj"^'^, j = 1,2 are bounded in L'^{U) independent of e. 
Therefore, extracting weakly convergent subsequences, not relabeled, we have that J^^^''^ •^j"'' Pj"''^ ^ 
P'"', and passing to the limit in (I4.22p and ()4.23|) we obtain 

div^f - VPf = /j- = divlff - VPf . 

and thus 

(6.10) div.Ff = V (Pf - Pf) + divT^/, J = 1, 2. 
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Similarly, ^f''^, PPI'^ are bounded in L'^{It x U) independent of e. Therefore, extracting weakly convergent 
subsequences as before, we obtain 



ipq 



(6.11) diyj^l"-^ = V (PJ'"'" - ) + divj; 
From (pHI)) . ((OT|) we deduce 

(6.12) div (.Fr'^e(ti;),,) = div iyi'e{w)j„) + V [(Pf'^ - Pf )e(ti;) 



Pr- - Pr Ve(«;),, + LfT'^ - . Ve(«,) 



and 
(6.13) 



iXiY {T^j"''" e{wt)pq) = div (^J^^'^e(u;f)pg 



V 



(Pf-~P7)e(«;0p, 



(Pf - Pf ) Ve(«;,)p, + (.^f - :^r) • ^<^t)p,^ 3 = 2, 3. 



Next, since div u,^ = 0, 
(6.14) 



(M^v 
(M^v 



(Pr^-Pf)e(«;; 



pq 



= 0, 



(Pr^-Pr)e(«;Op,J>=0, 

'^(^|'''-^'')(i-r)e(«;Op,(r)dr 



= 0. 



Integrating by parts and using (|6.12[) , (|6.13p and (|6.14p we have 



/(m'^,i(;'^) = -(w^div (^J^i''e(i(;)pg^ + div (^J^2'e(i(^t))9q)) 



(^i^div I y T)i {t -T + T)e{wt)pq{T)dTdxdt 



where 



^ -(«^ i^r^ - T\') ■ Ve{w)pq) - {{Tr'' - • Ve(«;0 



(^f'^-^r) (t-r + T). Ve(u;t)p9(T)dTda;dt) 
-(«^ (Pr^^ - Pf) . Ve(«;)p,) ~ (Pr-^ - pD • Ve(«;,)p,> 
^ (Pl^^^ - Pf) - r) • Ve(i«Opg(^)^^^^^a;^^i> 

since m*^ G L^{It, Hq{U)), we can apply Lemma [5731 which yields lim^^o = 0. Then we have (|6.9p . where 
the components of the effective tensors A{x), B{x), C{t, x) are defined by 



(6.15) 



A — T R — T C — T 



This completes the proof of the theorem 16. II 
Proof of the main theorem. 

To obtain (|7.12|) . we pass to the limit in (j2.13|) using Lemma [5.31 Next, insert — N'^w into (j2.14|) . 
The limit of the inertial terms is given in Proposition 15.41 and the limit of the term containing T'^ ■ e{w^) 
is provided by Theorem 16.11 Together, these results yield (|7.13p . The divergence-free constraint (|7.1ip is 
obtained by straightforward passing to the limit in div = 0. 
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7. Fluid-structure interaction 

Compared to the previous sections, the main difference now is lack of elhpticity in . In this section we 
assume ~ ^i^o- This means that phase one is a Kelvin- Voight viscoelastic material, and phase two is a 
Newtonian fluid. To deal with degeneration of we modify (|4.ip as follows. 

(7.1) N''w = w'' = w + J^ nP'^^^t-T + T)e{w)pq{T)dT + 

mP''^'{t^T + T)e{wt)pg{T)dT + W^'. 



Here, n^''^, m^^'^, 0^ are as in (|4.ip . (|4.3p . respectively. We note, however, that the auxiliary problems for 
^P9,e^^P9.e will be different. Differentiating (|7.ip we obtain 

(7.2) A'e (w') ~ B'e (wl) = ^r''e(ti;)p, + (.Ff^') (t^Opg + 



J Tl'^^'{t~T + T)e{w)„{T)dT + j Tr'\t-T + T)e{wt)p,{T)dT 

+A'{g\+gl)-B'{dtg\+dtg'^) 
+A'eiVc^')-B'eiV(bl), 



where 



J^P"-' = -B' {iP'i + e (m??'")) , 
TP^-^x) = A'e (nP«'=) - B'e (nf'") , 
.Ff '"(i, a;) = A'e (mf«'") - B^e (mf'") , 

(7.3) ^1 " ^ / (^""'^'(^ - T + T) ^e{w)pg{T) + {nPi-^t - r + T) ® Ve(-u;)p,(T))T) dr. 

(7.4) 3l = lJ^ {mP'^ 'it - T + T) ® Ve(«;t)pg(T) + (m?'«'^(t - r + T) ® Ve(ii;Opg(T))^) dr. 

In the statements of the following auxiliary problems we drop p, q to simplify notations. 

• First auxiliary problem. Find rij, € Hq{U) satisfying 

(7.5) div(74^Jf«-B"e(n|.))- VPf = div = 0, 

• Second auxiliary problem. Find G Hq(U) satisfying 

(7.6) -div(B''(7f«-f e(m|,)))- VP| = /2, div = 0, 

• Third auxiliary problem. Find n"^ G W satisfying 

(7.7) -div{{A' - B'dt)e{n'))-WPi ^ f^, div = 0, n^T) ^ n'j,. 

• Fourth auxiliary problem. Find G W satisfying 

(7.8) -div{{A^ - B'dt)e{m'')) -VPI = f^, div m'^ = 0, m^(r) = m|,. 

Since is still elliptic, the problems (|7.5p . (|7.6p can be analyzed exactly as problems in Section All 
the results in that section apply without change. The problems (|7.7p . (|7.8p were dealt with in Section HTT] 
The most important condition is still elhpticity of B^ , and the only change that is needed is in the proof of 
(i) in Proposition 14. 21 where we used elhpticity of A"^ . Now, to prove (i) note that 

(7.9) J J A'e{u)-e{u)dxdt+^ J B'e{u)-e{u)dx{t)= J {f ,u) H-i(^u),H^,{U){r)dT 
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holds for almost all t E It- Estimating the right hand side of (|7.9p we have 



(7-10) {f,u)H-i(u),Hl,{U)iT)dr < II / Wh-^u) W II u WH^i^u) W^^t 

< (^f II / Wl-HU) ir)dr^ {^[ II « Wlum i^)dr 



1/2 / „^ X 1/2 



<ll / \\l^{It,H-^U))\\ U \\l^{It.,H},{U)) 

Since the first term in the left hand side of (|7.9|) is non-negative, from (|7.9p . (|7.10p we deduce using ellipticity 
of 

^/3i sup / e(M) • e{u)dx{t) <\\ f \\lhIt,h-Hu))\\ « WL^ir.miu)) ■ 
This yields (i) (with a different constant i/JiT^^) after observing that 

W '^WhiiT m(U))= f f e{u) ■ e{u)dx < T sup [ e{u) ■ e{u)dx{t). 
Jo Ju teirJu 



All the arguments in Sections [51 El apply with minor changes due to the presence of four fluxes Tj, 
J = 1, . . . , 4, instead of three. The result can be summarized as a theorem. 

Theorem 7.1. In the case if fluid- structure interaction, the limits 'p,v, u satisfy 

(7.11) divtT=0, 

and the integral identities 

(7.12) / -pf^(j){(),x)dx- i i -pdt(l)dxdt~ I i pv ■ \/(j)dxdt = 0, 

Ju J It Ju J It Ju 

(7.13) — / PqVq ■ xpdx — / / pv ■ dtipdxdt — / / 'pv(^v-'S/ij}dxdt 

Ju J It Ju J It Ju 

+ T ■ e{xp)dxdt = 0. 

J It Ju 

for all smooth test functions <f), tp, such that div ■0 = 0, and (j),xp are equal to zero on dU and vanish for 
t>T. _ _ _ _ 

Moreover, there exist the effective tensors A G L'^{U),B G L^iU) and C G L'^{It x U),D e L'^{It x U) 
such that the effective deviatoric stress T satisfies 

(7.14) T = Ae{u) + Be{v)+ f C{t - T)e(u){T)dT + f D{t - T)e{v){T)dT 

Jo Jo 

Remark. The effective tensors are obtained as weak L^{It x U) limits of the four fluxes: 
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Appendix A. Moving interface in the inertial terms and frozen interface in the 

CONSTITUTIVE EQUATIONS 

In this section we present formal calculations leading to the weak formulation of the momentum balance 
equation. 

1. Inertial terms. 

(A.l) [ [ dt{p''v^)-'ipdxdt= [ [ 9'dtip'v')-ipdxdt^ 

J It Jv' J It Ju 

-I I 0%p''v^)-dtipdxdt- [ [ [p^v^ ■^|J)^te'dxdt- 
J It Ju J It Ju 

P^OqVq ■ ■0(0, x)dxdt. 



(A.2) / / dj{p'vlv'^)ii}dxdt= / / 9'dj{p''vlv'j)^^dxdt 

I It Jv^ J It Ju 



[ [ e^p^'vtv^jdjxp^dxdt- [ [ {p^v^ ■ip){v- ■\/0^)dxdt. 
J It Ju J It Ju 



IItJU JItJU 

Combining (jA.ip and (jA.2[) we obtain 

(A.3) ( [ [dt{p'v'')+dW{p^v^v)]-xpdxdt = 

J It Jv^ 

- I p^9oVa ■xp{0,x)dxdt- I I e^ip^v') ■ dti>dxdt- [ [ 9^ {p^v^) ■ dtipdxdt 

Ju J It Ju J It Ju 

[ [ {p'v' ■ip){dt9' ■Ve'')dxdt. 

J It Ju 

When 0^ satisfies the interface evolution equation, the last term in the right hand side is zero. If the interface 
were frozen, then this term would be present in the weak formulation of the momentum equation. Since 9^ 
is piecewise constant, the weak formulation would contain a non-physical term supported on the interface. 
To avoid such non-physical terms, one needs to use 0^ in the inertial terms of the momentum equation. 
2. Constitutive equation. Moving interface assumption combined with the Hook's law would lead to a 
non-physical dissipation of the elastic energy. Indeed, let the elastic part of the stress be written as 

(A.4) 9l{x)A^e{u') + [1 - 9l{x))A^e{u'), 

The important condition here is that 9q is independent of t. The stiffness tensors , A^ of the phases are 
supposed to be constant. Formally multiplying (jA.4p by v"^ and integrating by parts we obtain 

(A.5) / / {9oA^ + {l-9a)A^)e{u'')-e{dtu'')dxdt 

Jo Ju 

^ \ I f 9t [{OoA^ + (1 - 9o)A')e{u') ■ e{u')] dxdt 
^ Jo Ju 

(A.6) = i ^ {9oA' + (1 - 9o)A^)e{u^) ■ e{u^)dx \ ^ . 

This expresses the fact that the elastic energy changes by the amount of work done by elastic forces, with no 
dissipation. If one were to use 9' it, x) in (|A.5|) . differentiation in time would not commute with multiplication 
by 9^A^ + (1 - 9^)A^, and ((X6)) could not be obtained. 

Appendix B. Existence of weak solutions. Outline of the proof. 

In this Section we outline the proof of existence of global weak solutions for the system (|2.13p , (I2.14p for 
each fixed £ > 0. Since e is fixed, we drop superscript to simplify the notation. We follow closely pjj. Sect. 
2.3, 2.4. 

1. The initial conditions. The initial conditions satisfy (|2.15p . (|2.16p . 
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2. Formal a priori estimates l|2.20p - (|2.25p are obtained as explained in Sect. 12.21 In particular, renormal- 
ization as in [l^, Sect 2.3 is used to get \{x G U : a < p{t, x) < /3}\ for each < a < /3 < oo, where | • | 
denotes Lebesgue measure. This implies (|2.20|) . 

3. Compactness results. Since we wish to approximate exact solutions, a compactness result is needed. 
Suppose that we have two sequences satisfying the conditions < p" < C div ~ 0, a. e. on 
It X [/, II \\l^{It,hI{u))^ C, ^ V weakly in L'^{It, Hq{U)). Moreover, assume that 

9tp" + div(p"7;") = 

in V'{It X U), p"|t>"|2 is bounded in L°° {It , (U) , and we have 

for aU (f) e L«(/t, W"'''i{U)) such that div = 0. 

Finally suppose that the initial conditions for the density Pq satisfy 

Po Po 

in L^{U). 

Then, by the compactness Theorem 2.4 in TS] (it applies without any change), we have 
p" ^ pin C{[Q,T],LP{U)), for aU 1 < p < c5o, 

6p 



VP^r in L^ilT, L^'iU)), for 2 < p < oo, 1 < r < 



3p-4' 

V, in L\It, L^^U)), for 1 < 6* < 2 on the set {p > 0}, 

4. Construction of smooth approximate solutions. As in [18], Sect. 2.4, construct solutions {p,v) of the 
approximate system 

(B.l) dtp + div{pvs)^Q 

(B.2) dtipv) + div{pvs (g) v) - div {Ase{u) + Bse{v)) + VPs = in 

div t> = in P', 

where vs^ As^ B^, are smooth regularizations of the respective quantities. The initial conditions are regular- 
izations of the original ones. Then, using a fixed point argument as in Theorem 2.6, we can prove existence 
of smooth solutions to (jB.ip . (jB.2[) . Existence of a fixed point follows from the a priori estimates. The only 
issue that needs to be explained here is bootstrap regularity of the constructed solutions. The following 
proposition replaces Proposition 2.1 in [T5] 

Proposition B.l. Consider the system 

(B.3) cdtVi + h ■ Vvi — aAvi — mAui + diP = 0, 

div V — Q in It U , i — 1,2, i, with the initial conditions 

(B.4) t;(0,-) =-"0, w(0,-) = 0. 

Suppose that c £ L°°{It x U), a,m £ L°°(/t, W^^'°°(C/)), b £ {It , L°° {U)) , c> k, a> k, m> k a.e. on 
It 'x U for some k > 0; Vq £ Hq{U). Also, assume that a,m are independent oft. 

Then the system HB.S]) . has a unique solution {v,P) such thatv £ L'^{It, Il'^{U))nC{[0,T], II^{U), 

dtv £ L'^{It X U), VP £ I^{It X U). 

Outline of the proof. Compared to the proof of Proposition 2.1 in [18j . we have one new term mAui. 
Multiplying this term by dtVi and integrating by parts we have 

mAuidtVi ~ dt I Vui-mVvi— I m|Vwi|'^ + / Vui-VrndtVi 
u Ju Ju Ju 

Multiplying (jB.Sp by dfVi using the above identity, integrating by parts and summing over i we find 
[ c\dtv\^ + ^dt [ a\Vv\'^ + dt [ Vu ■ mVv [ m\Vv\^ ^ [ f, 

Ju 2 J{j Ju Ju Ju 
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where f = b ■ VvdfV — Vit • (Vto (g) dtv). Integrating this identity with respect to t we obtain 

nt 



\dtv\' + kl 
m|VM||V'u|(i) 



\S/v\'it) 



n\Vv\'+ f I f + \ I a\Wvm. 
Jo Ju ^ Ju 

j^m\Vu\\Vv\{t) <|1 m U^.(u) \Vv\\t) + ^ \\/u\\t) 



Next we write 



where we choose v — k/2. The term containing / is handled similarly, putting v in front of Jq Ju \dtv\ 
Combining the previous two inequalities with the standard a priori bounds on Vit, Vt> we have 



where Ck depends only on k and the data. This yields a priori estimates on v in L°°{It, Hq{U)) and on dtV 
in L^{It X U). 

Now we can write (IB. 31) as 



(B.5) A{av + mu) - VP = h, div v ^ 0, in U, 

V g Hq{U),u e Hq{U) for almost all t € /t. Also, h is bounded in L^{It x J7) in terms of the data. Next, 
estimating pressure P exactly as in [T^, Prop. 2.1, we conclude that av + mu £ L'^{It, H'^{U)). Since a, m 
are smooth and positive, this implies 



(B.6) 



dtu H u = g, 



where g G L^{It, H^{U)). Now, formally, u{t, •) = e " '^^g{T, ■)dT, and hence 



{t)< 



d'g {r)dT ] <t d'a {T)dT < T / 8' a {T)dT. 



Integrating over U we deduce that u is bounded in L°°{It, H'^{U)) in terms of the data. Then using 
V = —^u + g we obtain that v is bounded in L^{It, H'^{U)) in terms of the data. The formal calculations 
can be easily justified by an approximation argument. 

■ 

5. Passage to the limit. This is done using compactness from step 3 exactly as in [TS] . 
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